Charmonium dynamics in Au+Au collisions at y/s = 200 GeV 
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The formation and suppression dynamics of J , Xc and ^' mesons is studied within the HSD 
transport approach for Au + Au reactions at the top RIfIC energy of — 200 GeV. Two prominent 
models, which have been discussed for more than a decade, are incorporated, i.e. the 'hadronic 
comover absorption and reformation' model as well as the 'QGP threshold' scenario, and compared 
to available experimental data. Our studies demonstrate that both scenarios - compatible with 
experimental observation at SPS energies - fail severely at RHIC energies. This combined analysis 
- together with the underestimation of charm elliptic flow - proves that the dynamics of c, c quarks 
are dominated by partonic interactions in the strong QGP (sQGP) and can neither be modeled by 
'hadronic' interactions nor described appropriately by color screening alone. 
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According to current understanding, the evolution 
of the universe in the 'Big Bang' scenario has pro- 
ceeded from a quark-gluon plasma (QGP) to color neu- 
tral hadronic states within the first second of its life- 
time. In this context, the dynamics of ultra-relativistic 
nucleus-nucleus collisions at Super-Proton-Synchrotron 
(SPS) and Rclativistic-Heavy-Ion-GoUider (RHIG) ener- 
gies are of fundamental importance as reflecting the prop- 
erties of hadronic/partonic systems at high energy den- 
sities. The c, c quark degrees of freedom are of particular 
interest with respect to a phase transition from baryonic 
matter to the QGP, since cc meson states might no longer 
be formed in the very hot fireball due to color screen- 
ing [U, S| • This initial intuitive expectation has guided 
experimental studies for almost two decades. However, 
more recent lattice QGD (IQGD) calculations have shown 
that the J/'i> survives up to at least 1.5 Tc (Tc « 170 to 
185 MeV) such that the lowest cc state remains bound 
up to rather high energy density [1, d, Q . On the other 
hand the Xc and appear to melt soon above Tc. It 
is presently not clear, if also the D and D* mesons will 
survive at temperatures T > but strong correlations 
between a light quark (antiquark) and a charm antiquark 
(quark) are likely to persist [7]. One may speculate that 
similar correlations survive also in the light quark sec- 
tor above Tc, such that 'hadronic comovers' - most likely 
with different spectral functions - might show up also at 
energy densities above 1 GeV/fm'^, which is taken as a 
characteristic scale for the critical energy density. 
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The production of charmonium in heavy- ion collisions, 
i.e. of cc pairs, occurs dominantly at the initial stage 
of the reaction in primary nucleon-nucleon collisions. At 
the very early stage the cc pairs are expected to form 
color dipole states which experience i) absorption by in- 
teractions with further nucleons of the colliding nuclei 
(c/. Refs. These cc color dipolcs can be absorbed 

in a 'pre-resonance state' before the final hidden charm 
mesons or charmonia (J/^, Xc, are formed. This 
absorption ~ denoted by 'normal nuclear suppression' - 
is also present in p + A reactions and determined by a 
dissociation cross section as ~ 4 to 7 mb. Those char- 
monia or 'pre-resonance' states that survive normal nu- 
clear suppression during the short overlap phase of the 
Lorentz contracted nuclei furthermore suffer from ii) a 
possible dissociation in the deconfined medium at suffi- 
ciently high energy density and iii) the interactions with 
secondary hadrons (comovers) formed in a later stage of 
the nucleus-nucleus collision. 

In the QGP 'threshold scenario', e.g the geometrical 
Glauber model of Blaizot et al. [ll| as well as the perco- 
lation model of Satz (sj, the QGP suppression 'ii)' sets 
in rather abruptly as soon as the energy density exceeds 
a threshold value Sc, which is a free parameter. This 
is motivated by the idea that the charmonium dissoci- 
ation rate is drastically larger in a quark-gluon-plasma 
(QGP) than in a hadronic medium On the other 
hand, the extra suppression of charmonia in the high 
density phase of nucleus-nucleus collisions at SPS ener- 
gies [l2, has been attributed to inelastic comover scat- 
tering (c/. H, M, m, H [13, M, [H and Refs. therein) 
assuming that the corresponding J/^'-hadron cross sec- 
tions are in the order of a few mb [20i, i22, 23]. In 
these models, 'comovers' should not be viewed as asymp- 
totic hadronic states in vacuum but rather as hadronic 
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FIG. 1: The J/^ nuclear modification factor Raa (El for Au+Au collisions at ^/s = 200 GeV as a function of the number of 
participants Npart in comparison to the data from lOj for midrapidity (full circles) and forward rapidity (full triangles). HSD 
results for the 'QGP threshold melting' scenarios are displayed in terms of the lower (green solid) lines for midrapidity J/'i/'s 
{\y\ 5: 0.35) and in terms of the upper (orange dashed) lines for forward rapidity (1.2 < y < 2.2) within different recombination 
scenarios (see text). The error bars on the theoretical results indicate the statistical uncertainty due to the finite number of 
events in the HSD calculations. Predictions for the ratio i3;_i^(^'')a^//_B^^ (>^/*I')c"j/* as a function of the number of participants 
Npart for Au+Au at ^/s = 200 GeV are shown in the lower set of plots. 



correlators (essentially of vector meson type) that might 
well survive at energy densities above 1 GeV/fm^. Ad- 
ditionally, alternative absorption mechanisms might play 
a role, such as gluon scattering on color dipole states as 
suggested in Refs. [13, H^, 1131 or charmonium dissoci- 
ation in the strong color fields of overlapping strings [1^ . 

We recall that apart from absorption or dissociation 
channels for charmonia also recombination channels such 
as D + D ^ Xc + meson {Xc = (J/^',Xc,*') play a 
role in the hadronic phase. A previous analysis within 
the HSD transport approach [23, [s^l - employing the co- 
mover absorption model - demonstrated that the char- 
monium production from open charm and anticharm 
mesons indeed becomes essential in central Au+Au col- 
lisions at RHIC. This is in accordance with independent 
studies in Refs. [13, [1^ and also with the data from 
PHENIX [U. On the other hand, the backward chan- 
nels - relative to charmonium dissociation with comoving 
mesons - {Xc+ meson D + D) were found to be prac- 
tically negligible at the SPS energies. 

In the present study we extend our previous investi- 
gation [34I within the 'comover model' and the 'QGP 



threshold scenario' to the energy of -^s = 200 GeV and 
compare to the PHENIX data. The questions we aim 
at solving is: 1) can any of the models be ruled out by 
the present data sets and 2) do the recent PHENIX data 
provide a hint to a different dynamics of charm quarks 
at top RHIC energies? 

The explicit treatment of initial cc production by pri- 
mary nucleon-nucleon collisions is the same as in Ref. |32l | 
(see Fig. 1 of Ref. [H] for the relevant cross sections) and 
the implementation of the comover model - involving a 
single matrix element Mq fixed by the data at SPS ener- 
gies - as well as the QGP threshold scenario are as in [s^ . 
Consequently no free parameters enter our studies below. 
We recall that the 'threshold scenario' for charmonium 
dissociation is implemented as follows: whenever the lo- 
cal energy density e[x) is above a threshold value e^, 
where the index j stands for J/^,Xc,^', the charmo- 
nium is fully dissociated to c + c. The default threshold 
energy densities adopted are ei = 16 GeV/fm'^ for J/^E", 
£2 = 2 GeV/fm^ for Xc and £3 = 2 GeV/fm^ for and 
provide a fair reproduction of the data at SPS energies 
(except for in the 'threshold scenario'). The reader is 
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FIG. 2: Same as Fig. [T] for the 'comover absorption scenario' 
including the charmonium reformation channels without cut 
in the energy density (l.h.s.) and with a cut in the energy 
density Scut = 1 GeV/fm'^ (see text for details). 

referred to reference [13] for details. 

The energy density £{r;t) - which is identified with 
the matrix element T'^^{r;t) of the energy momentum 
tensor in the local rest frame at space-time (r, t) - be- 
comes very high in a central Au-|-Au collision at ^/s — 200 
GeV, according to the HSD calculations, where baryons 
with approximately projectile or target rapidity are omit- 
ted. In the center of the reaction volume, e{r;t) initially 
reaches values well above 30 GeV/im^ and drops below 
1 GcV/fm"^ roughly within 5-7 fm/c. We recall that in 
HSD explicit hadronic states are allowed to be formed 
only for e{r;t) < 1 GeV/fm^. 

In the theoretical approach, we calculate the J/^* sur- 
vival probability ^j/* and the nuclear modification fac- 
tor Raa as 
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dN{J/-^)AA/dy 
NcoU-dN{J/^)pp/dy' 



(1) 
(2) 



where N^l^ and denote the final number of J/^ 

mesons and the number of J/^'s produced initially by 
BB reactions, respectively. Note that Nj[^ includes the 
decays from the final Xc- In (0), dN{J/'^)AA/dy denotes 
the final yield of J/^E* in AA collisions, dN{J /'^)pp/dy is 
the yield in elementary pp reactions, while Ncoii is the 
number of binary collisions. 



We start with a comparison of Raa{JI'^) © for 
Au + Au collisions as a function of the number of partic- 
ipants Npart to the data from in the upper part of 
Fig. [T] The results for the 'threshold melting' scenario 
(without the reformation channels D+D — + ( J/^, Xc, 
-|- meson) are displayed on the l.h.s. of Fig. 1 in terms 
of the lower (green) solid line for midrapidity J/'^'s 
{\y\ < 0.35) and in terms of the upper (orange) dashed 
line at forward rapidity (1.2 <|y| < 2.2). The exper- 
imental data from PHENIX |lO| are given by the full 
circles at midrapidity and by triangles at forward rapid- 
ity. In this simple scenario, practically all charmonia are 
dissolved for Npart > 50, due to the high energy densi- 
ties reached in the overlap zone of the collision, which 
is clearly not compatible with the PHENIX data and 
indicates that charmonium reformation channels are im- 
portant. Here we explore two scenarios for charmonium 
reformation: a) we adopt the notion that hadronic cor- 
relators (with the quantum number of hadronic states) 
survive above Tc and the reformation and dissociation 
channels {D + D ^ (J/'^^ Xc ^'') + meson) are switched 
on after a formation time Tf = 0.5 fm/c (in the local rest 
frame) and b) the hadronic states are assumed to persist 
only below s{r;t) < 1 GeV/fm^ and thus the reforma- 
tion and dissociation channels {D + D ^ (-^/^i Xc, + 
meson) are switched on only for energy densities below 
1 GeV/fm^. The results for the model a) are displayed 
in the upper middle part of Fig. 1 and demonstrate that 
for Npart > 200 an approximate equilibrium between the 
reformation and dissociation channels is achieved. How- 
ever, here the calculations for forward rapidity match the 
data at midrapidity and vice versa showing that the ra- 
pidity dependence is fully wrong. Furthermore, the J/^* 
suppression at more peripheral reactions is severely over- 
estimated. The results for the model b) are shown in 
the upper right part of Fig. 1 and demonstrate that the 
dissociation and reformation channels no longer reach an 
equilibrium even for most central collisions. The J/^ 
suppression as a function of centrality as well as rapidity 
is fully off. Summarizing our model studies, we have to 
conclude that the 'threshold melting + reformation sce- 
nario' is incompatible with the PHENIX data and has to 
be ruled out at top RHIC energies. 

In the lower parts of Fig. 1, we show the results for the 
ratio of the and J/^* dilepton yields (given by their 
cross sections multiplied by the corresponding branching 
ratios) which have no experimental counterpart. Here 
the two recombination models give finite ratios as a func- 
tion of centrality but predict a larger to J/vP ratio 
at forward rapidity than at midrapidity which is a con- 
sequence of the higher comover density at midrapidity. 
Experimental data on this ratio should provide further 
independent information. 

The ratio RAAiJ/^) in the 'comover -I- recombination 
model' is displayed in the upper part of Fig. [2] in compar- 
ison to the data from using the same assignment of 
the lines as in Fig. 1. The l.h.s. shows the results for the 
'default' comover reformation and dissociation channels 
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(as in Ref. [s^]) whereas the r.h.s. corresponds to the re- 
sults when the comover channels are switched on only for 
energy densities e{r;t) < Scut — 1 GeV/fm'^. The latter 
scenario shows a suppression pattern which is in strong 
contrast to the data both as a function of Npart and ra- 
pidity. The default scenario (l.h.s.) gives a continuous 
decrease of Raa{J /^) with centrality, however, an oppo- 
site dependence on rapidity y due to the higher comover 
density at midrapidity. The to J/^* ratio is displayed 
in the lower parts of Fig. 2 and shows a decreasing ra- 
tio with centrality similar to the results at SPS energies 
[3^ . However, independent from experimental results on 
this ratio, the 'comover -|- recombination model' has to 
be ruled out at RHIC energies, too. 

In concluding and summarizing our study, we have in- 
vestigated the formation and suppression dynamics of 
J/^*, Xc and 'f' mesons - within the HSD transport ap- 
proach - for Au-f An reactions at top RHIC energies of 
^/s — 200 GeV. Two controversial models - discussed 
in the community for more than a decade, - i.e. the 
'hadronic comover absorption and reformation' model as 
well as the 'QGP threshold melting scenario', have been 
compared to the available experimental data from the 
PHENIX Collaboration When adopting the same 

parameters for cross sections (matrix elements) or thresh- 
old energies as at SPS energies [s^l, we find that both 



scenarios - compatible with experimental observation at 
SPS energies - fail severely at RHIC energies and can 
safely be excluded. This provides a clear answer to the 
question 1) raised in the introduction. 

We point out, furthermore, that the failure of the 
'hadronic comover absorption' model goes in line with 
its underestimation of the collective flow V2 as well as the 
underestimation of Raa{pt) of leptons from open charm 
decay as investigated in Ref. [331 • This strongly sug- 
gests that the dynamics of c, c quarks are dominated by 
partonic interactions in the strong QGP (sQGP) which 
can neither be modeled by 'hadronic' interactions nor de- 
scribed appropriately by color screening alone. This also 
gives an answer to question 2) of the introduction. 

Since the open charm suppression is also underesti- 
mated severely in perturbative QCD approaches, the na- 
ture of the sQGP and its transport properties remain an 
open question (and challenge). 
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